Creating collective many-body states with highly excited atoms 
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We study the collective excitation of a gas of highly excited atoms confined to a large spacing 
ring lattice, where the ground and the excited states are coupled resonantly via a laser field. Our 
attention is focused on the regime where the interaction between the highly excited atoms is very 
weak in comparison to the Rabi frequency of the laser. We demonstrate that in this case the 
many-body excitations of the system can be expressed in terms of free spinless fermions. The 
complex many-particle states arising in this regime are characterized and their properties, e.g. their 
correlation functions, are studied. In addition we investigate how one can actually experimentally 
access some of these many-particle states by a temporal variation of the laser parameters. 

PACS numbers: 67.85.-d, 32.80.Ee, 42.50.Dv, 71.45.-d 



I. INTRODUCTION 

Ultracold atoms provide a unique toolbox to study 
many-particle physics under very clean and well-defined 
conditions. The precise control over their interactions 
and their trapping potentials allows to study the dy- 
namics of phase transitions as well as the preparation 
of strongly correlated quantum states [l[ . 

While so far the majority of experiments is carried out 
with ground state atoms, exploiting the unique proper- 
ties of highly excited states is gradually moving into the 
focus of experimental and theoretical efforts. Atoms in 
highly excited states can interact strongly, i.e., the in- 
teraction strength can be of the order of several tens of 
MHz at a distance of several micrometers. The corre- 
sponding quantum dynamics takes place on a microsec- 
ond timescale and thus is orders of magnitude faster than 
the atoms' external dynamics. Such scenario is usually 
referred to as 'frozen gas' @, A number of experi- 
mental groups have studied the excitation dynamics of 
such system using Rydberg states of alkali metal atoms 
[3| which were excited from an ultracold gas. Here a dra- 
matic reduction of the fraction of excited atoms was ob- 
served once the atomic density was too high or the inter- 
action between excited states was too strong [j| @ . This 
is a manifestation of the so-called 'Rydberg blockade' Q] 
effect that is responsible for the collective character [8[ 
of Rydberg excitations in dense gases. 

Very recently the power of Rydberg states to estab- 
lish a controlled interaction of single atoms trapped in 
distant traps has been demonstrated in a series of im- 
pressive experiments H M, ED, El • Strongly supported 
by these results, highly excited atoms nowadays are be- 
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lieved to have a manifold of applications ranging far be- 
yond traditional atomic physics. Indeed, exploiting the 
properties of atoms in Rydberg states permits the study 
of spin systems at criticality [1 31 ] , the quantum simulation 
of complex spin models [14[ , the investigation of the ther- 
malization of strongly interacting many-particle systems 
[IH and also the implementation of quantum information 
protocols 

In a recent work (Ref . [13] ) we showed that the unique 
properties of Rydberg atoms allow the creation of en- 
tangled many-particle states on a one-dimensional ring 
lattice on a short time scale. Finding simple ways for 
creating entangled many-particle states is of importance, 
since such states have a number of applications, e.g, they 
serve as resource for the creation of single-photon light 
sources [l8|. for improving precision quantum measure- 
ments [19] and for measurement based quantum informa- 
tion processing. 

In this paper we will go into depth and largely ex- 
pand on our previous study. We show that excited many- 
particle states of a laser-driven gas of Rydberg atoms on 
a ring lattice can be obtained analytically in the limit of 
strong laser driving. We give a detailed derivation of the 
system's Hamiltonian in Sec. HU The construction of the 
many-body excitations, their eigenenergies and their cor- 
relation properties are analyzed in Sec. Mil In Sec. llVlwe 
discuss thoroughly how these states can be excited in an 
experiment. We conclude with a summary and outlook 
in Sec. M 



II. THE SYSTEM 

We study a gas of bosonic ground-state atoms confined 
to a deep large spacing optical or magnetic [13, HH ring 
lattice with periodicity a « /im (see Fig. [T|). The Wan- 
nier functions \I/fc(x) are localized at the k-th site with a 
width (7 <C a. We assume the external dynamics of the 
atoms to be frozen, i.e., no hopping and hence no par- 



2 



A" l R >* 
Q 

— IPX 



k-1 

c 

k 



o _ k+1 

C C c 



FIG. 1: Ring lattice with spacing a being much larger than 
the extension a of the Wannier functions (deep lattice). The 
internal atomic degrees of freedom at each site are described 
by the (collective) states \P) k and \R) k , coupled by Q. 



tide exchange between the lattice sites is present. This 
is well justified as the internal (electronic) dynamics - 
in which we are interested here - takes place on a much 
shorter timescale, of the order of hundred nanoseconds. 
We consider two electronic levels which are denoted by 
\g) and \r). Here \r) is a Rydberg ns-state which - due 
to its quantum defect - is well isolated from any other 
electronic level. It is coupled to the ground state \g) via 
a laser with Rabi frequency f^o and detuning A. Within 
the rotating wave approximation, the Hamiltonian de- 
scribing the coupling of the atoms to the laser field reads 
(with h = 1) 



fc=i 



rik, 



(1) 



k=l 



where b 



(b k and r k ) represent the creation (anni- 



hilation) of a ground and a Rydberg state, respectively, 
and n k = r k r k stands for the number of atoms in state |r) 
at the k-th site. We will consider throughout this paper 
the case where each lattice site is occupied by the same 
number of atoms, Nq. This is achieved, for example, if 
the system is initialized in a Mott-insulator state. 

The interaction between the Rydberg atoms is given by 
the van-der-Waals potential Cq/x 6 , that is quickly decay- 
ing with the distance x between atoms. Nevertheless, as 
Cq scales with the eleventh power of the principal quan- 
tum number n, the interaction can strongly affect the 
excitation dynamics of atoms that are separated by sev- 
eral micrometers. This strong interaction gives rise to the 
so-called blockade effect 0, HH ■ We consider a scenario in 
which the simultaneous excitation of two or more atoms 
to the Rydberg state on a single lattice site is blockaded. 
Thus, on each lattice site k, only the two states 



\P) k = 

\R) k = 
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are accessible, where S is the symmetrization operator. 
The effective Rabi frequency for the laser coupling be- 
tween these so-called (super)atom states (see Fig. [1]) is 



given by = flo^/No. Taking all this into account, in Eq. 



(P) we can replace Slo r k b k 
M * Jk) 



n 



, A - (n/2)[ai k) +ia { y % 

where a x ^' and ay" 1 are the Pauli spin matrices. 

Since a <C a (see Fig. [T]) we can rewrite the van-der- 
Waals potential between two (super)atoms in the state 
\R) located d sites apart as (3d = C^jx\, where Xd is 
the separation between those sites. As already pointed 
out, (3d is quickly decaying with the distance. In par- 
ticular, the next-nearest neighbor interaction is a factor 
of 64 cos 6 (7r/L) smaller than the nearest neighbor one 
(02 = /V(2cos(tt/£)) 6 ). We will thus only focus on 
the nearest neighbor interaction which is well-justified 
for large enough lattices. The interaction Hamiltonian 
for the entire atomic ensemble, with f3 = (3\ 1 then reads 



L 
k=l 



with the Rydberg number operator n k = [1 + ai k ^]/2 and 

the boundary condition = ■ 

In summary, the complete Hamiltonian that drives the 
dynamics of our system can be written as 



Eh 



+ An k + f3n k n k+ i 



(2) 



fc=i 



The system can be described as a periodic arrangement 
of spin-1/2 particles, where the two spin states, corre- 
sponding to the two internal states of the (super) atoms, 
\P)u and \R) k , interact via an Ising-type potential. In 
this picture, the Rabi frequency f2 and the combination 
of A + (3 can be effectively interpreted as perpendicular 
magnetic fields. Hence, the relevant parameters in our 
system will be: a) the ones related to the laser, i.e., the 
single-atom Rabi frequency ilo and detuning A, which 
can be time-dependent and b) the interaction between 
Rydberg atoms represented by (3. 



A. Constrained dynamics 

Throughout this paper, we consider the regime where 
the detuning is much smaller than both the collec- 
tive Rabi frequency (laser driving) and the interaction 
strength, i.e., |A| <C ^,/3- As a consequence, the be- 
havior of the system will be determined by the ratio of 
the latter two parameters. Here we focus on the limit 
f2 ^> (3, i.e., the laser coupling is much stronger than 
the interaction between atoms. In this regime the first 
term of the Hamiltonian is the dominant one and it 
is convenient to make it diagonal by means of a rotation 
of the basis. This is achieved by the unitary transfor- 
mation U = Ilfc=i CX P ( — *f CT y fc ') which brings a x — * a z 
and a z — > —o x . When applied to our Hamiltonian @, it 
yields 



H 



U^H Bpin U 



(3L 



+ H XY + Hi + H 2 , 



(3) 



3 



with 
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H 2 = 
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where H xy is the famous icy-model of a spin chain with 
a transverse magnetic held. 

Let us now analyze the importance of the individual 
contributions of H . As we can see in Fig. [21 the spectrum 
of H decays into manifolds of states which are separated 
by gaps whose width is approximately 2Q. This is caused 

by the dominant first term of H xy , i.e., ^X^fe 17 -^- The 



eigenstates of a. 
- given by 



(fc) 



are - in terms of the (super) atom states 



\±)k = ^[\P) k ±\R) h 



with ai k) |±k = ± |±) 



Thus, each of the manifolds 



that determine the coarse structure of the spectrum is 
spanned by a set of product states that have the same 
number of (super)atoms in the state |+). In Fig. [2J these 
manifolds are denoted by m, which is the eigenvalue of 
the states with respect to the operator ai k \ 

The second term of H xy conserves the total number 
of |+) (super)atoms. In other words, it couples only 
states that belong to the same m-manifold and that are 
nearly degenerate. As a consequence, the strength of 
these intra-manifold couplings due to H xy is proportional 
to (3. Conversely, Hi and H 2 couple states that belong 
to manifolds with different number of (super)atoms in 
the state |+). In particular, Hi and the last term of H 2 
flip one of the (super) atoms from |+) to |— ) or viceversa. 
Thus, the coupled states belong to different manifolds 
with Am = ±1, energetically separated by 2f2. The two 
first terms of H 2 drive a similar process, flipping always 
two contiguous (super) atoms in the same state simul- 
taneously, i.e., |++) — * | ) or | ) — * |++). As a 

result, these terms connect states with eigenvalue m to 
those with m ± 2, which are separated roughly by 4f2. 
These features are reflected in Fig. [2] The transition 
rates between m-manifolds corresponding to Hi and H 2 
can be estimated by second order perturbation theory to 
be of the order A 2 /fi and (3 2 /£l, respectively. Hence, for 
sufficiently strong driving CI 3> /3, their contribution can 
be neglected and the system's dynamics is constrained 
to the m-manifolds. As a consequence, the Hamiltonian 
that drives the intra-manifold dynamics, H xy , effectively 
drives the dynamics of the entire system in this parameter 
regime. This Hamiltonian is analytically solvable, and we 
thus have access to the actual spectrum and eigenstates 
of the system. The diagonalization of this Hamiltonian 
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FIG. 2: Level structure in the regime Q S> /3 and |A| <C n. 
The spectrum splits into manifolds which can be labeled by 
the quantum number m of the operator ^ fc ai k ^ . For suffi- 
ciently large Q, the coupling between manifolds that is estab- 
lished only by Hi and H2 can be neglected. The (constrained) 
dynamics inside the m-subspaces is then determined by H xy . 



relies on the so-called Jordan- Wigner transformation and 
a Fourier transform that we explain thoroughly in the fol- 
lowing paragraph [23| . 



B. Jordan- Wigner transformation on a ring 

The Pauli matrices in the Hamiltonian ^ obey anti- 
commutation and commutation relations when they be- 
long to the same and different sites, respectively. Thus, 
the algebra is neither bosonic nor fermionic. This dif- 
ficulty can be overcome by the Jordan- Wigner transfor- 
mation, 



fe-i 

n 



k-l 



Cfc = 



n H°) - 



(*) 



(7) 



which introduces the operators ct and Ck that obey the 
canonical fermionic algebra 

{c|, cj} = S itj {c|, c]} = {a, c 3 } = 0. 

After this transformation, the Hamiltonian ([4]) takes on 
the form 



L r 

E 

k=l 



2Q 



ctcfe 



[c\ci + c\c L 



(4 



1 



Cfc+i + c[ +1 c k 



(8) 



Thus, the Hamiltonian has been transformed into one 
which describes a chain of spinless fermions with near- 
est neighbor hopping. The last term of Hamiltonian ([S]) 
appears due to the periodic boundary conditions. It de- 



pends on the operator n_ 



- T L etc 



which counts the 



total number of fermions, which is also equivalent to the 
number of (super)atoms in the state |+). Thus, depend- 
ing on the parity of the number of fermions of the state, 



4 



-ff X y reads 



the diagonalized Hamiltonian ([9]) reads 



L 

e • 

fc=i 



= > '2fi | c< k c k 



!A Z^Wt t 

2 +4L (4 c fc+i + 4+i c fc 



/c=i 



T^ ( C L C 1 + c \ c l) , 

for even (e) or odd (o) parity, respectively. 

These two cases can be accounted for simultaneously in 
a convenient way by introducing a matrix representation 
for the fermionic operators. They are projected onto the 
subspaces with even and odd eigenvalue of n + by means 
of the projectors P c/o = [l ± e i7m +] /2, with P c + P a = 
1. Since the Hamiltonian H xy conserves the number of 
fermions, i.e., [i? xy , e Z7r ™+] = 0, it is diagonal in this 
representation and can be decomposed as 



-LQ 



E 



2Q 







AlA r , 



(10) 



As we will see in the next section, the introduction of the 
matrix- valued fermionic operators has the advantage that 
excited states can be constructed by applying products 
of AJj to the ground state. As a consequence, this matrix 
notation allows us to automatically distinguish between 
the odd and even fermion cases, which otherwise has to 
be done manually. 



III. MANY-BODY STATES 
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We now introduce new matrix-valued creation and anni- 
hilation operators of the form 



ll 



4 
4 



7fe 



c k 
c k 



which obey the fermionic algebra provided c k and c\. are 
fermionic operators. The Hamiltonian can now be con- 
veniently written as 



with 



20 



L 

E 



fe=i 



L-l 



k=l 



(7I71 + IiIl) e 



(9) 



1 

-1 



The diagonalization of the Hamiltonian © is achieved 
by performing the following Fourier transform 



7fe = — 2^ VnkA ™ ^ k 

n=l 

with the Fourier coefficients 

Vn. ' 



v n— 1 





■ 2w 1 



e 

The operators A^ n and A„ are matrix-valued 



A„ = 



VnPo 
VnPc 



with 77^ and r\ n being fermionic creation and annihilation 
operators, respectively. Defining the eigenvalue matrix 

e n ss 



e„ = 2 



[^(n-1/2)] 




COS |_T7 n ] 



A. Symmetries 

The symmetry properties of our system impose cer- 
tain selection rules for the excitation of the many-particle 
states. In order to understand this, let us start our 
analysis of the excited states by studying the symme- 
tries of the Hamiltonian. Because of the special arrange- 
ment of the sites, the Hamiltonian ([2]) and also (fTOf are 
invariant under cyclic shifts and reversal of the lattice 
sites. This can be formally seen by representing these 
two symmetries through the operators X and TZ, respec- 
tively. Their action on the spin ladder operators are 

*V£ fc) * = 4 fe+1) and TZ^a^TZ = a^' k+1 \ from where 
follows that [H spin , X) = [-ff spin , TZ] = 0, i.e., both of them 
correspond to conserved quantities. Thus, if the system 
is initialized in an eigenstate with respect to X and TZ, 
the time evolution will not take place in the entire Hilbcrt 
space, but merely in the subspace spanned by the states 
with the same quantum number with respect to X and 
TZ. 

This observation is highly relevant for our system. In 
practise, the natural initial situation will be that in which 
all atoms are in the ground state, i.e., |0) = Ilfe=i \P)k- 
This state has the above-mentioned properties, i.e., it is 
invariant under cyclic shifts and the reversal of the sites: 
Aft |0> = |0) and TZ^ |0) = |0). We will refer to such a 
state that has eigenvalue 1 with respect to X and TZ as 
being fully-symmetric. Hence, only the states from this 
fully-symmetric subspace can be actually accessed in the 
course of the system's time evolution under Hamiltonian 
@. In the following we will thus focus on constructing 
excited states that belong to this subset. 



B. Fully-symmetric states 

The ground state of Hamiltonian (jTUJ) is given by 



i G >=n 



k=l 



5 



and it is fully-symmetric. Excited states that contain 
TV fermions are in general formed by successive applica- 
tion of the creation operator, i.e., \N pq ...) — A|Aj . . . \G). 
However, not all combinations will give rise to states that 
belong to the fully-symmetric subset. 

Let us start considering the possible cases of a single- 
fermion excitation. For a fully-symmetric state we re- 
quire Ot |l p ) - Ot A t | G ) = t A t \G) = A| \G) =\l p ), 
i.e., 

OUtO \G) = A| |G> , 

with O being a placeholder for X and 1Z. After some 
algebra one finds that 



Since e ivn + \G) = \G) and e i7rc i Cl \G) = \G), only the 
single excitation with p ~ L is symmetric under cyclic 
shifts and reversal. Hence, the only one-fermion state 
that can be reached by the time-evolution reads 

|1)=A[|G). 

To have a better physical understanding of this state, it 
is convenient to write it in terms of the atomic operators, 



ID 



i L 



\G). 



Thus, |1) is a spin wave or, in other words, a superatom 
that extends over the entire lattice. These states are of 
interest since they can be used as a resource for single 
photon generation. 

For the two-fermion states, we follow the same proce- 
dure and demand 

Ot A t A tO|G)=AtAt|G). 

One finds that 

T^t]\t]\TI = r — 



VL-q+lVL-p+l 



+ 



1) 



From this, one sees that the condition p+q— 1 = L has to 
be accomplished. As a result, the fully-symmetric states 
are 



K4.- P+ i \G) 



with p = 1 . . . [L /2J . These are entangled states formed 
by superpositions of two-atom excitations in the ring with 
opposite momentum. This is more clearly seen by writing 
everything in terms of the Pauli matrices 



sm 



k>k' 



2n 



b-l/2)(fc-fc')Vf^ fe,) l G > 



These states are potentially interesting for the production 
of photon pairs. How they can be actually accessed will 
be discussed in Sec. IIVI 

Finally, let us illustrate how the three-fermion excita- 
tions arc formed. We have 



'lL- 



-e i2 ^ p+q+r H-A-A-r^ n+ - 



and thus fully symmetric three-fermion states are of the 
form 

|JW> = i= (At AjAt - A{_ p A{ _ q A{ _.) \G) , (11) 

with p + q + r = L, 2L. Writing these eigenexcitations 
back in terms of the spin operators yields 

\3pqr) = ~J~3/2 X/ X/ £ P1 r 
k>k'>k" P crm 
(pgr) 



2tt 



(kp + k'q + k"r) 



a^a^ap \G) 



where s pqr is the Levi-Civita symbol. In a similar way, 
states with higher number of fermions are obtained. 



C. Energy spectrum 

Now that we have analyzed the eigenstates of the sys- 
tem we will focus on the corresponding eigenenergies. In 
the course of this investigation we will also perform a 
comparison of the analytic results to the ones obtained 
from a numerical diagonalization of the Hamiltonian ([2]). 
This will allow us to assess the accuracy of our analytical 
approach. 

Let us begin with the ground state energy. From Eq. 
(fTUl) we can read off the value 



E (: = -l [ n-| 



(12) 



where we have included the general energy-offset f3L/4 
(see Eq. ©). For A = 0, Q = 10, j3 = 1 and L = 
10, the result is Eq = —97.5. This is to be compared 
with the numerical value of —97.63 which is obtained by 
diagonalizing the Hamiltonian |2]). We find both results 
to be in good agreement. For the first excited state we 
obtain 

E l = E G + 2fl+^. 

Using the same set of parameters, the energy of the 
single-fermion state is E\ = —77.0, which is very close 
to the numerically exact value —77.11. The energies of 
higher eigenexcitations are given by 



E 2p = E G + 4fi + /3cos 



6 




Q=10 



FIG. 3: Spectrum of Hamiltonian H Q for a lattice of L = 10 
sites versus the laser driving f2 in units of f3. In the right 
insets, the energies of the two- and three-fermion states are 
shown for Q. = 10. Five two-fermion and eight three-fermion 
eigenenergies arise as it is analytically predicted for this lattice 
size. 
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-58.09 


-58.17 
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-58.45 


-58.54 



10, A = 0, Q = 10 and (3=1 and comparison with the 
numerically exact values. 



with p = 1 . . . [L/2\ , for the two-fermion case and 



E 



3pqr 



£ G + 60 + ^ 



2tt 
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with p + q + r = L, 2L, for the three-fermion one. For 
L = 10, we obtain five and eight different eigenenergies 
for the two- and three-fermion states, respectively (see 
insets in Fig. [3]) . In the Tables U and |TT] we perform 
a comparison between the analytical and the numerical 
results. A difference of less than a 0.2% is observed in all 
cases. 



The discrepancies between analytical and numerical 
values are mainly caused by second order energy shifts 
due to H\ and H% (Eqs. ([5|) and ([6])). These contribu- 
tions vanish only in the limits /3 /r2 — > and A/O — * 0. 
Here, we will calculate them for a finite ratio. There is a 
constant term in Hi which is proportional to A that gives 
rise to a global energy shift = LA/2. Being aware of 
this shift facilitates the comparison between the numer- 
ically exact and the approximate analytical eigenvalues 
for A f 0. 

Let us focus first on the ground state. Hi and H2 
only couple states whose number of fcrmions differ by 
one or two (Fig. [2]). As a consequence, only the states 
|1) and 

the energy of the ground state. It yields 



2 P ) contribute to the second order correction of 



E. 



(2) 



L\A + /3\' 



80 + 2/3 

LL/2J 

E 

P =i 



4 



sin 2 [^(p-1/2)] 



40 + (3 cos 



1/2)] 



Analogously, we calculate the energy shift of the first 
excited state, |1), due to Hi and H^. In this case, we 
have to compute the effect of the states \G), \2„) and 
|3„ 



Jpqr/ 



r(2) 



The resulting energy correction is given by 
L|A + /3| 2 I 



80 + 2/3 

LL/2J 

< E — 
U 20 



L 

cot 2 \l(p- 



1/2)] 



2 cos [f (p-l/2)]-l) 



For the parameters A = 0, O = 10, j3 = 1 and L = 10, 
these shifts yield e£ } = -0.14 and e{ 2) = -0.10. The 
corrected energies of the ground and the single-fermion 



state are now Eq = E ( 



(0) 



Ei = E{ 0) 



e[ 1] + EY> = -77.10, much closer to the 
numerically exact ones of —97.63 and —77.11, respec- 
tively. We will later see that these energy corrections 
can be useful for the selective excitation of many-particle 
states in the lattice. 
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D. Correlation functions 
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Numerical 
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10 


-36.19 
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-36.75 
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-37.10 


-37.15 
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-37.36 
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6 


-37.65 


-37.71 
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10 


-37.81 


-37.87 
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-38.00 


-38.05 
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3 


5 


-38.00 


-38.06 



TABLE II: Energies of the eight three-fermion states |3 p «j r ) 
for L = 10, A = 0, f2 = 10 and (3 — 1 and comparison with 
the numerically exact values. 



In this subsection we are going to study the density- 
density correlation function of the many-particle states. 
This quantity measures the conditional probability of 
finding two simultaneously excited atoms at a distance 
x from each other normalized to the probability of un- 
corrected excitation. It is defined - for a fully-symmetric 
state |\&) - as 

ff2{x,V) = a 1, 

where we have used (n a )^, = (n;,)^ for all sites. The 
correlation function will give gi(x, "J") = when two sites 



7 



separated by a distance x are completely uncorrelated, 
and g2(x, > (< 0) for correlation (anticorrelation) 
between the sites. 

In particular, for the case = \2 p ), g2(x,2 p ) can be 
analytically calculated. In terms of the expectation val- 
ues of the spin operators, the correlation function reads 

92{x,2 p ) = (<t<V 1+x) +a W a^ +x) ) 2 . For ^Owe 
have 32(0, 2 P ) = 1 and for x > the calculation yields 



9i{x, 2 p 



4 

L2 



(i — 2a; ) cos 



2tt 



(p - l/2)a: 



-2 sin 



T (p- 1/2)1 



cot 



?(p-l/2) 



By inspecting this expression for the allowed values p = 
1, . . . L^/2J , some general statements can be made: 

a) Independently of the total number of sites L, there 
are always two 'extremal' cases (see Fig. |4^) which cor- 
respond to p = 1 and p = [L/2\: for p = 1, the cor- 
relation function shows a positive maximum at x = 1, 
i.e., nearest neighbor, and then decreases monotonically 
and smoothly with the distance, staying always positive; 
for p = [L/2\, the nearest neighbor is pronouncedly an- 
ticorrelated, the next-nearest neighbor is correlated and 
this pattern of correlation-anticorrelation persists with 
increasing distance. 

b) For p = [L/2\ , the oscillations of #2(2;, 2 P ) are more 
pronounced for L than for odd L, see Fig. 2^. The 
ratio of the amplitudes of the correlations for x — 1 and 
x = [L/2\ is, 



92(1,2$) 

<?2(f,2 # ) 



1 



.92(1,2^1) 



92 ( 



£-1 
2 > 



J 



in the even and odd cases, respectively. Also, for an even 
number of sites, the correlation functions of the two ex- 
treme cases accomplish g2{x,2h) = (—l) x g2(x,2\), i.e., 
the envelope of the oscillating function g2(x 1 2i L ) is given 
by the smoothly decreasing g%(x, 2\). 

c) For a fixed value of p, the amplitude of the corre- 
lations decreases with increasing number of sites as 1/L, 
as can be seen in Fig. 0b. 

Numerically, we have observed agreement to the ana- 
lytical results shown in Fig. [5] As expected, this agree- 
ment improves with a decreasing ratio /3/f2. The corre- 
lations could be directly monitored experimentally pro- 
vided that a site-resolved detection of atoms in the |+)- 
state is possible. The next section will deal with the open 
question of how these correlated states can be experimen- 
tally accessed. 



IV. EXCITATION OF MANY PARTICLE 
STATES 

Our aim is to selectively excite correlated many-body 
states by a temporal variation of the laser parame- 
ters. Initially the atoms shall be in the product state 




T*T-»-d. 










^Js^b. 2;-°--© 




| f I -s-L=20 
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1 ~L=24 
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FIG. 4: Density-density correlation functions of the \2 P ) 
states, a: For p = 1 and p — [L/2\, the correlations show 
completely different behavior, i.e., smoothly decreasing and 
strongly oscillating, respectively. These oscillations are much 
more pronounced for the even value of L = 24 than for the 
odd, L = 25. b: The magnitude of the correlations decreases 
as the number of sites L is enhanced, as can be seen for 
L = 20,24, 28. 



|0) = Ilfc=i \P)k ano - t ne laser shall be turned off, i.e., 
f2o(0) = and A(0) = Ao- Starting from these initial 
conditions, the goal is to vary flo(t) and A(i) such that 
at the end of the sequence, i.e., at t = £fi na ,i, the detun- 
ing is zero and the laser driving is much larger than the 
interaction (A(tfi na i) = and 0(tfi n al)//3 3* 1). This fi- 
nal situation corresponds to the right-hand side of the 
spectrum presented in Fig. [3l 

Once a desired many-particle state has been populated, 
and due to the limited lifetime of the highly excited lev- 
els which is in the order of several fis (e.g., 66 (is for Rb 
in the 60s state), we want to map it to an stable con- 
figuration. To do so, we first turn off the laser (tt = 0) 
and then switch on a second one whose action can be 
described by the Hamiltonian 



L 

E 

k=l 



■0 slnkrik+i. 
fe=i 



(13) 



In this expression, s k and Sk stand for the creation and 
annihilation operators of an single-atom stable storage 
state \s) on site k, respectively. In the limit where the in- 
teraction is much smaller than the Rabi frequency of this 
transition, i.e., Q s 3> /3, we can neglect the second term 
of this Hamiltonian. Thus, performing a global 7r-pulse 
to the considered many-particle state means to perform 
the mapping \r) — > |s), such that a stable configuration 
is achieved. 

Hence, the difficulty lies in finding a 'trajectory' or 
sequence (Q,o(t), A(t)) for which at t = £fi na ) only a single 
many-particle state is occupied. We propose two different 
methods in the following. 



8 




100 



FIG. 5: Possible trajectories (A(t),f2o(i)) through the spec- 
trum of H with A / (units of (3). a: When A < 0, the 
ground state at Q, = does not coincide with the initial state, 
|0), and the energy of the initial state goes through a number 
of avoided crossings. A possible path through them to reach 
the state |523456) is shown, b: If Ao > 0, the initial state |0) 
is adiabatically connected to the ground state \G). 



A. Direct trajectory 

In certain cases, one can guess a trajectory 
(A(t), Oo (t) ) like the ones shown in Fig. [5] that eventu- 
ally connects |0) with a desired eigenstate of H xy (24l.[25|. 
but this is not always possible. The general appearance 
of the laser sequence strongly depends on the sign of the 
initial detuning A . In Fig. [5] the two possible scenar- 
ios (taking A ^ 0) are depicted. For A < 0, the initial 
state is not the ground state of the system when the laser 
is turned off (f2 = 0). As a consequence, this initial state 
suffers several avoided crossings with other levels when f2 
is increased. Thus, it is not easy to find a path through 
the spectrum that connects it to a single desired eigen- 
state of -ffxy, as the one shown in Fig. [3^. A more general 
framework for finding a proper trajectory is provided by 
Optimal Control theory [2a |. Here, the desired fidelity 
with which the final state is achieved can be set and cer- 
tain constraints on the trajectory can be imposed. This 
method is successfully applied to quantum information 
processing 27], molecular state preparation [28[ and op- 
timization of number squeezing of an atomic gas confined 
to a double well potential [291 ]. The case of A > will 
be treated in the next subsection. 



B. Excitation from the ground state 

We present in this work a different route to populate 
single many-particle states. This is accomplished in two 
steps: First, one has to prepare the ground state \G) of 
Hamiltonian ([J| in the limit £1 ^> /?; once the ground 
state is populated, the single-fermion and two-fermion 
many-particle states can be accessed by means of an os- 
cillating detuning, that gives rise to a time-dependent 
Hi. 

Step 1: Let us start by explaining how to vary the 



laser parameters to prepare the ground state \G) . In 
particular, when setting A > 0, the ground state of 
the system at 17 = coincides with the initial state |0). 
With increasing SI, it is adiabatically connected to the 
ground state \G) of H xy (see Fig. [5Jd) . The problem that 
we can encounter here is that non-adiabatic transitions 
to other energy levels occur when increasing $7, so that 
we do not populate only \G) but also other states. To 
avoid this, we choose a large enough value of A when 
the laser is still turned off (f2 = 0). This increases the 
energy gap between |0) and other energy levels, and, as a 
consequence, suppresses non-adiabatic transitions. This 
initial detuning can be decreased as Q increases so that 
in the desired regime, i.e., fldnai = ^(ifinal) /3, it is set 
to zero. As an example, we propose the following shapes 
of n(t) and A(t) 



n(t) 

A(t) 



^finai sin 

A :i 



/ TTt 



sin 



2 ( T< 



inal 



(14) 
(15) 



The obtained fidelity | (3>(ifi n ai)|G) | for different values 
of the initial detuning Ao and time intervals ifmai is given 
in Fig. [SI where ^(ifinai) stands for the wavefunction of 
the final state. It is actually possible to populate the 
desired state with high fidelity, e.g., over 99% is achieved 
for all considered lattice sizes with Ao = 45/3 and ifi na i = 
0.9/3 -1 . We find that: i) the fidelity depends only weakly 
on the lattice size although the dimension of the Hilbcrt 
space grows exponentially with L, and ii) as expected, for 
a fixed value of the initial detuning, the fidelity increases 
with the increasing length of the time interval. Note 
that the timescale of this whole process is limited by the 
lifetime of the Rydberg state. 

If there is only one atom per site, and based on the 
fact that |G) = rifc=i l — )& i s a product state, an alter- 
native procedure to this adiabatic passage can be envis- 
aged. Starting from the vacuum |0) (also a product state 
with every atom in \g}) 7 we perform a global 7r/2-pulse 
to the single-atom transition \g) — > \s). As a result, we 
obtain a product state where every atom is in a super- 
position [\g] + i \s}] j\[2. In a second step, the 7r-pulse 
with the mapping laser described by the Hamiltonian 
(fT3]) and with fl s ^> f3, transfers every atom to the state 
[\g) — \r}] /V%, i-e., we have prepared the ground state 
\G). It is worth remarking that this method eliminates 
the lifetime limitation in this first stage. 

Step 2: Let us show now how to address the single- 
fermion and two-fermion states from this ground state 
\G). As we explained in section lTl A[ the Hamiltonian Hi, 
associated with the detuning, drives transitions between 
neighboring manifolds, i.e. Am = ±1, (see Fig. [2]). We 
exploit this fact and introduce an oscillating detuning of 
the form A(t) = A osc cos (tut) . If we tune lu to coincide 
with the gap between two given states, this detuning acts 
effectively as a laser that couples them resonantly with a 
Rabi frequency that is proportional to A osc . 
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FIG. 6: Fidelity | ($(tfi na i) \G) | 2 when populating the ground 
state of H xy from the initial state via variation of the param- 
eters of the laser Q(t) and A(t) in the form given by (|14p and 
(|15|) . respectively. Several initial values of the detuning and 
time intervals, as well as different lattice sizes, are considered. 
For a fixed value of Ao (units of /3), better fidelities are ob- 
tained for larger time intervals (units of 1//3). For a fixed 
time interval, there is an optimal value of Ao for each size of 
the lattice, around Ao ~ 45. 



Using this oscillating detuning, we want to transfer 
the population from the ground to the first excited state 
(Fig. Uh)- To do so, io is tuned to be on resonance with 
the corresponding energy gap, i.e., u = u>i = E\ — Eq, 
and by a 7r-pulse we populate |1). One has to take into 
account that in the limit of Q, ^> /3 the energy gap be- 
tween any two neighboring manifolds is equal, i.e., also 
higher lying excitations are populated. To avoid this ef- 
fect and address only the |1) state, we can choose a not 
too large value of Q. In this regime, the second order level 
shifts caused by Hi and H2, that are roughly given by 
A 2 /fi and f3 2 /Q, respectively (see section ITlI CI) , become 
increasingly important. In particular, as it is sketched in 
Fig. [7^, the gap between |1) and any of the \2 p ) levels 
becomes more and more different from uj\ and, as a con- 
sequence, the unwanted transitions fall out of resonance. 
Analogously, the same procedure could be used to ad- 
dress the two-fermion many-particle states (see Fig. [JJd). 
The first 7r-pulse resonant with the \G) — > |1) transition, 
is followed by another 7r-pulse with u> tuned to coincide 
with the energy gap of the specific |1) — > \2 P ) transition, 
uj = ll>2 = Ei v — E\ . The separation between neighboring 
|2 p ) states is of the order of f3 and the Rabi frequency of 
the transition is proportional to A osc . As a consequence, 
to populate only a single level of the two-fermion mani- 
fold, the parameters have to accomplish that (3 ^S> A osc 
and, at the same time, A osc has to be large enough in 
order to perform the transfer at a time interval that is 
much shorter than the lifetime of the Rydberg state. 



V. CONCLUSIONS AND OUTLOOK 

In this work we have studied the collective excitation 
of a laser-driven Rydberg gas confined to a ring lattice. 



FIG. 7: Sketch of the excitation of the single-fermion and two- 
fermion states by means of an oscillating (radiofrequency) de- 
tuning using a not too large value of fi. a: In a first step, the 
population is transferred by a 7r-pulse to the single-fermion 
state by tuning the frequency of the detuning on resonance 
with the gap ui = u)\. b: A second 7r-pulse with u tuned to 
match LU2 addresses the corresponding 2 P ) state, bearing in 
mind that /3 S> A osc in this step. 



We have focused on the regime in which the interaction 
between the highly excited states is much weaker than 
the laser field. We found that the corresponding system 
can be described as a chain of spinless fermions whose dy- 
namics is driven by the xy-moAe\. This Hamiltonian can 
be analytically solved and, by exploiting the symmetries 
of the system, we were able to completely characterize 
the many-particle states arising. In particular, we have 
shown that the first excited state of the Hamiltonian cor- 
responds to a spin wave or to an excitation which is com- 
pletely delocalized all over the lattice. The two-fermion 
states could be expressed as a superposition of excitation 
pairs and an investigation of their density-density cor- 
relation function has been performed. We have demon- 
strated that the qualitative behavior of these correlations 
differs substantially from one state to another of the same 
two-fermion manifold, going from a smoothly decaying 
function to a pronounced correlation-anticorrelation pat- 
tern. The analytical eigencnergics of the xy-Hamiltonian 
were compared to the numerical exact ones of the com- 
plete Hamiltonian, and excellent agreement between both 
results has been found. Finally, we have investigated 
several paths for the selective excitation of the many- 
particle states. One of them relies on the variation of 
the laser parameters with time, finding trajectories from 
the initial to a given final many-body state. The other 
possibility we have presented makes use of an oscillating 
detuning which allows to access excitations starting from 
the ground state of the Hamiltonian. In each step, a tt- 
pulse is performed with the frequency of the oscillation 
matching the energy gap between the involved states. 

Throughout this work we have considered an homoge- 
neous occupation of the sites of the ring lattice. The situ- 
ation of having a randomly fluctuating number of atoms 
per site would effectively lead to a disorder potential for 
the fermions, as outlined in Ref. [l7|. This would im- 
ply as well a change in the symmetry properties of the 
system, so that more states become accessible by a time- 
evolution (e.g., L possible single-fermion states instead 
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of only the fully-symmetric one). In addition, we have 
assumed that the atoms are strongly localized, a -C a 
(Fig. [I]). Taking into consideration the finite width of 
the wave packet would lead to another kind of disorder, 
this time associated to the interaction parameter [3. 

As we have pointed out, the main problem one has 
to face in this system is the limited lifetime of the Ry- 
dberg states, which is in the order of several microsec- 
onds. One could think of preparing a parallel system to 
the one described in this work but using polar molecules 
[30l l3lj , to overcome this lifetime limitation. An inter- 
esting extension is also the investigation of the system 
in two-dimensional geometries, e.g., triangular or square 



lattices, as well as several rings disposed in concentric or 
cylindric configurations. In all these cases, the symme- 
tries of the particular arrangement of the sites might give 
rise to new interesting many-particle states. 
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